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ABSTRACT 

We argue that gravitational microlensing is a feasible technique for measuring the 
mass function of brown dwarf stars in distant galaxies. Microlensing surveys of the 
bulge of M31, and of M87 in the Virgo cluster, may provide enough events to differ- 
entiate the behaviour of the mass function of lenses below the hydrogen burning limit 
(though we find that M87 is a more favourable target). Such objects may provide a 
significant supply of baryonic dark matter, an interesting possibility for the study of 
galactic dynamics. Furthermore, these systems have different metallicities than the 
solar neighbourhood, which may affect the mass function. These considerations are 
relevant in the context of star formation studies. 

Key words: gravitational lensing — stars: low mass, brown dwarfs — 
cosmology:observations — dark matter — galaxies:haloes 



1 INTRODUCTION 

Gravitational microlensing, by which dark objects are de- 
tected by their magnification of bright sources being mon- 
itored, is a rapidly growing field, with many applications 
in astrophysics. Using microlensing to find dark objects in 
the halo of the Milky Way was first proposed by Paczynski 
(1986), and elaborated by Griest (1991). Several groups 
have extended the possibilities of microlensing by the so- 
called pixel technique (Crotts 1992; Baillon et al. 1993), 
where source stars are not resolved due to the high degree 
of crowding. However, successive images can be analysed, 
usually by taking image differences, to uncover the bright- 
enings due to variability in the individual sources, some of 
which will be due to microlensing. This technique has been 
applied successfully by several collaborations (Tomaney & 
Crotts 1996; Alcock et al. 1999a,b; Wozniak 2000). Previ- 
ously, it has been difficult to extract physical information 
from these pixel microlensing lightcurves because of a de- 
generacy in the lightcurve shape with magnification. The 
method of Gondolo (1999) for measuring the optical depth 
in pixel microlensing without knowledge of the magnification 
in an event has been extended to measurements of the mass 
function of lenses (Baltz 2000). In this paper we show how 
to apply this technique to measurements of the low-mass 
end of the stellar mass function, as has been done in classi- 
cal microlensing surveys for brown dwarfs (Zhao, Spergel & 
Rich 1995) and planetary- mass objects (Alcock et al. 1998) 
in our own galaxy. We hope to extend such limits to distant 
galaxies as well. 



2 BROWN DWARFS 

Brown dwarfs are defined to lie in the mass range between 
the hydrogen-burning limit (75 Mjup) and the deuterium 
burning threshold (13 Mjupi IMq = 1047Mjup). A theoret- 
ical review of these objects is given by Chabrier and Baraffe 
(2000). Radial velocity searches have detected both brown 
dwarfs and giant planets from 75Mjup to 0.25Mjup in or- 
bits around solar-type stars. The observational situation for 
brown dwarfs is reviewed by Basri (2000). Studies of young 
star-forming regions have recently found evidence for iso- 
lated objects in the giant planet mass range (5 — 15Mjup) 
(Zapateo-Ossorio et al. 2000). Isolated brown dwarfs have 
also been reported in various star-forming regions (Lucas & 
Roche 2000; Hillenbrand & Carpenter 2000; Najita, Tiedo & 
Carr 2000). The frequency of such objects is not known, but 
they seem to be sufficiently numerous that the IMF must 
certainly continue from the hydrogen burning limit to the 
giant planet regime, and could either be flat, exhibit a weak 
power law decline, or even rise below the brown dwarf limit. 

Observations of isolated brown dwarfs and giant planets 
are necessarily focussed on searching in star- forming regions, 
via the technique of infrared imaging. We show here that 
gravitational microlensing provides a potentially powerful 
technique for exploring the initial mass function of brown 
dwarfs in distant galaxies. 

We investigate several possible mass functions for this 
type of object, below the hydrogen burning limit. Our pro- 
posed method will not be able to distinguish fine features in 
the mass function, so we investigate several toy mass func- 
tions. It is the gross features of the mass function that we are 
immediately concerned with, and simple model mass func- 
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tions will suffice. We consider the class of mass functions 
(taking M in units of Mq, and with arbitrary normalisa- 
tion) 



log 



dN 



dlog M 



[A log A/ + 5(1- log M)2] 



(1) 



where A and B are free parameters (note that 5 = cor- 
responds to a pure power-law mass function dN/d log M oc 
M-^). These mass functions agree at ten solar masses, and 
taking A = 1.35, have the Salpeter (1955) slope at ten solar 
masses. This prescription ensures that the mass functions 
produce equal light, as most of the visible light from stellar 
populations comes from massive stars, heavier than the sun. 
The uncertainty in the population thus lies in the mass to 
light ratio, not in the total light, which is as it should be. 

Taking B = 0.25 gives a mass function with approx- 
imately the Salpeter (1955) slope down to O.O2M0. This 
model has a large amount of stellar mass below the hydrogen 
burning limit, indicating a large amount of baryonic dark 
matter locked up in brown dwarf stars. Taking B = 0.32 
gives a mass function turning over at around O.OSMq, simi- 
lar to the Miller-Scalo (1979) mass function. Lastly, we take 
B = 0.48, giving a mass function turning over at around 
OAMq, similar to that of Gould, BahcaU & Flynn (1996). 
This model has the least amount of baryonic dark matter. 




Figure 1. Self-lensing event rate towards M87 bulge with ACS 
on HST. The sample spacing is four times daily. There is a much 
lower sensitivity to lenses below about 10~^ A^q. We illustrate 
delta function mass functions, together with the B = 0.32 Miller- 
Scalo type mass function. 



3 SURVEYS OF DISTANT GALAXIES 

Pixel microlensing surveys of large distant galaxies are quite 
attractive for learning about both the dynamical properties 
of the systems and the populations of objects from which 
they are made. We will discuss two possible targets to il- 
lustrate the capabilities of this technique with ground and 
space based telescopes. 

A pixel microlensing survey of M31, the Andromeda 
Galaxy, is quite feasible using a three-meter class ground 
based telescope. This large spiral galaxy is about 725 kpc 
distant, and is the nearest large galaxy to the Milky Way. 
We are primarily interested in the bulge of this galaxy, as it 
affords the densest star fields and the largest optical depth 
for star-star lensing. However, finite source size effects for 
this target would require nearly continuous monitoring for 
reasonable sensitivity at low masses to be achieved, so we 
will not consider M31 further. 

Using a space-based telescope, the reach of the pixel 
technique is much longer. The giant elliptical galaxy M87, at 
the center of the Virgo cluster at a distance of about 15 Mpc, 
is well within the reach of a pixel microlensing survey by the 
Hubble Space Telescope (HST) . We take a mass model based 
on the work of Tsai (1993), 



p(r) = 3.8 



1 + 



kpc 



M, 



PC 



1, 1 + 0.275 log ( 

Vkpc 



(2) 



(3) 



The Advanced Camera for Surveys (ACS) is an ideal 
instrument for this type of observation. The future NGST, 
planned as an eight meter class space telescope, will allow 
a much more thorough survey, with better statistics. For 



both telescopes we assume exposures of thirty minutes (for 
the HST ACS, this is possible in one orbit), taken every six 
hours. More frequent sampling is punishing in that fewer 
events are detected with the same telescope resources (say 
30 orbits for HST), though smaller masses can be probed. 
Less frequent sampling may allow more events to be detected 
with the same telescope resources, but at the cost of the low- 
mass sensitivity that we desire for this program. For our 
purposes, six-hour sampling is about ideal, though with the 
NGST there is more room to take a larger sample spacing. 
We show the event rates for these two telescopes below, in 
Figs. I and | respectively (Baltz & Silk 2000). Rates for 
the NGST assume a field of view 4' square, with 50% more 
throughput than the HST ACS, and nine times the collection 
area. In all cases in this paper, we assume that seven samples 
above 2a must be collected, as Criteria A of Alcock et al. 
(2000). 



4 MEASURING THE MASS FUNCTION 

For each pixel microlensing event, we make in effect two mea- 
surements. These are the flux increase at maximum AFmax, 
and the full- width at half maximum time, tfwhm- As Wozniak 
& Paczyhski (1997) have clearly shown, even in the classical 
microlensing case, the shape of the lightcurve does not give 
much more information than these two parameters. How- 
ever, in the classical microlensing case, an additional mea- 
surement, that of the unlensed flux of the source star, is also 
made, allowing an accurate determination of the magnifica- 
tion of the event. 

The Einstein times of events have been used to esti- 
mate the optical depth to microlensing, and also to estimate 
the masses of the lenses. However, to make a good estimate 
of the Einstein time from a microlensing event, both ffwhm 
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Figure 2. Self-lensing event rate towards M87 bulge with NGST. 
The sample spacing is four times daily. The sensitivity drops con- 
siderably below about 10""^ A/q. The mass functions are the same 
as in Fig. ^ 



and the magnification are required. We thus seek a differ- 
ent characteristic of events. The time-scale tfwhm is prob- 
lematic, as it strongly depends on the unknown magnifica- 
tion. We instead follow Gondolo (1999) and use the quantity 
ujf = A_Finaxifwhm , which is effectively the Einstein time 
multiplied by the flux of the source star. This quantity is 
easily measured from pixel microlensing events. We choose 
to further normalize by dividing by the surface brightness 
fluctuation flux [F <x lO""'^™), 



^F^aayi ^fwhm 



F 



(4) 



The SBF flux F is measured by characterizing the pixel-to- 
pixel variations of the surface brightness of a galaxy image 
(Tonry & Schneider 1988). This normalisation gives the ap- 
proximate relation, with (3 being the event's minimum im- 
pact parameter in Einstein units, 



tp = tplQ 



-0.4(m — m) 



(5) 



V i 72 

when /3 ^ 1, or equivalently, at high magniflcation. Gon- 
dolo (1999) showed how to measure the optical depth to 
microlensing using this time-scale. We will use an extension 
of that method, described in more detail in Baltz (2000), to 
constrain the mass function of the lenses, in this case brown 
dwarfs. The method is based on producing a weighted his- 
togram of event rate with flux time-scale tp. Differentially, 
we want to study the dimensionless rate 

dT 



N, 



tp 



d log tp ' 



(6) 



given a delta-function mass function. The quantity Np then 
encodes the response of a microlensing survey to lenses of 
a given mass, effectively the number of events expected per 
decade in time-scale when monitoring for a time equal to 
the event time-scale. Equivalently, this quantity is closely 



related to the optical depth contributed by lenses producing 
events at a specific time-scale. 

For comparison, we can form quantities analogous to 
Np using the other relevant time-scales, namely the Ein- 
stein time (Ne) and the fuU-width at half maximum time 
(A'fwhm)- Their definitions exactly mimic Eq. ^ 

At fixed mass, one would naively think that the Einstein 
time is most useful, as the distribution is the narrowest of 
the three. However, the Einstein time can be measured only 
rarely. Furthermore, we find that the shape of the rate distri- 
bution with the Einstein time is sensitive to the exact details 
of the cuts used to define events. Interestingly, the shape of 
the rate distribution with the flux time is much less sensitive 
to the exact definition of an event, which is clearly a desir- 
able feature. This happens because changing the definition 
of an event usually amounts to changing the minimum value 
of the flux increase at maximum. This will simply change the 
short time-scale end of the tp distribution, leaving the peak 
shape intact in most cases. This point is further explored in 
Baltz (2000). Thus, it seems that the time-scale tp is per- 
haps preferable to the Einstein time in extracting the param- 
eters of the source-lens system. In Figs. |^ and ^, we display 
the function Np, given delta- function mass functions. Effec- 
tively, these are the smoothing kernels over which the true 
mass function can be measured. In Fig. ^ we plot the func- 
tion Np, based on the Einstein time-scale, for the same exact 
parameters as for the ACS in Fig. ^ As is clearly evident, 
the function Np is more sensitive near its peak than Np to 
changing the effective threshold, as is the case for changing 
the sampling strategy. Thus the discussion in this paragraph 
is validated: even if tp, and thus Np, were measurable, Np 
is preferred in cases where there is a high sensitivity to the 
threshold for events. 

The universality of the shape is broken by two effects: fi- 
nite source size and finite time sampling. The primary effect 
of the finite source size is that the magnification is bounded 
by the fact that the source stars are not point sources. This 
effect is obviously more pronounced for smaller Einstein 
rings, and thus for smaller masses. A finite time sampling 
implies a minimum event time-scale that can be detected. 
Lower mass lenses of course produce shorter events, thus 
more of the low-mass events are missed. The finite total ob- 
serving time would also subtract very long events, but for 
our purposes, this is not a concern, as we are interested in 
the short events due to 0.01-0.1 Mq stars. 



5 DISCUSSION 

We now illustrate the full Np distribution, integrated over 
the simple mass function discussed in Sec. 2. As discussed 
previously, Np for a single mass acts as a smoothing filter 
for the mass function. We show the results for the HST ACS 
in Fig. I and the NGST in Fig. ^ The pure Salpeter mass 
function has a peak at the shortest time-scale, the flattened 
models (Miller-Scalo and Salpeter with cutoff) has a longer 
peak, and the GBF-type model has the longest peak time- 
scale. With these surveys, the statistics should be sufficient 
to distinguish these models. We should state that the nor- 
malisation alone is probably insufficient to distinguish mod- 
els, but the position in time-scale of the peak of the rate 
distribution is quite robust. In fact, we have adjusted the 
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Figure 3. Dimensionless rate Np for ACS on HST. The curves 
are, left to right, for lenses of a single mass: 1O~*'~^'~^'~^'*''^M0. 
Samples are taken every 1.5 hours in the top curves, with the 
lower curves sampling every six hours. With six hour sampling, 
the sensitivity begins to decline at masses less than W~^Mq, and 
is quite low below 10~^Mq. 
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Figure 4. Dimensionless rate Np for NGST. The curves are anal- 
ogous to those in Fig. ^. With six hour sampling, the sensitivity 
only begins to seriously decline at masses less than 10~^Mq. 



normalisation of the B = 0.3 and B — 0.75 models so that 
the peak rates agree, in order to compare the time-scales. 
For the surveys we consider, the peak time-scale differs by 
a factor of 2-3 between a mass function that is flat for sub- 
stellar masses (Miller-Scalo), and one that that is sharply 
decUning at those masses (GBF). We note that mass func- 
tions that differ only below 0.1 Mq will be quite difficult 
to distinguish, but for mass functions which begin to differ 



Figure 5. Dimensionless rate Np for ACS on HST. The curves 
are analogous to those in Fig. ^, and in fact represent the same 
microlensing survey, with the time-scale tp rather then tp. 

at around a solar mass, the microlensing technique is quite 
powerful. 

We have shown that pixel microlensing can be a pow- 
erful tool for measuring the mass function of low mass and 
brown dwarf stars, less massive than the sun. Since this tech- 
nique is effective to very large distances, we have a chance 
to learn something about the universal properties of brown 
dwarf mass functions. 
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Figure 6. Mass functions with ACS on HST. We illustrate the 
measured Np for the three basic mass functions. The normalisa- 
tions of the B = 0.25 and B = 0.48 curves have been adjusted to 
compare the peak positions with the B = 0.32 case. 
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Figure 7. Mass functions with NGST. We illustrate the measured 
Np for the three basic mass functions. The normalisations have 
been adjusted as in Fig.H 
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